26/05/0^ 



Large spin corrections in TV^ = 4 SYM sl(2): still a linear integral equation 



i 



Diego Bombardelli, Davide Fioravanti and Marco Ross 

Sezione INFN di Bologna, Dipartimento di Fisica, Universitd di Bologna, 

Via Irnerio 46, Bologna, Italy 

Abstract 

Anomalous dimension and higher conserved charges in the s/(2) sector of = 4 
SYM for generic spin s and twist L are described by using a novel kind of non- 
linear integral equation (NLIE) . The latter can be derived under typical situations 
of the SYM sectors, i.e. when the scattering need not depend on the difference 
of the rapidities and these, in their turn, may also lie on a bounded range. Here 
the non-linear (finite range) integral terms, appearing in the NLIE and in the 
dimension formula, go to zero as s — > oo. Therefore they can be neglected at 
least up to the 0{s^) order, thus implying a linear integral equation (LIE) and a 
linear dimension/charge formula respectively, likewise the 'thermodynamic' (i.e. 
infinite spin) case. Importantly, these non-linear terms go faster than any inverse 
logarithm power (Ins)"", n > 0, thus extending the linearity validity. 
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1 Introduction 

The AdS/CFT duality [1] conjectures the equivalence between a string theory on the 
curved space-time AdSs x in the strong coupling and a conformal quantum field theory 
on the boundary of AdSs in the weak coupling. In particular, type JIB superstring 
theory should be dual to A/" = 4 Super Yang-Mills theory (SYM) in four dimensions. 
In particular, it would relate energies of string states to anomalous dimensions of local 
gauge invariant operators of the quantum field theory. In this context, the discovery of 
integrability in both free string theory and planar field theory was a great achievement, 
both from the conceptual and the practical (i.e. calculative) point of view, being the 't 
Hooft coupling, A = Sn'^g'^, the only non-running constant on stage (the string tension 
is proportional to g). Actually, integrable models appear as spin chain like Bethe 
equations, satisfied by 'rapidities' which parametrise on the one side the quantum string 
states (and their energies) and on the other side the corresponding composite operators 
(and their anomalous dimensions) in SYM, respectively. Actually, the initial result, 
which re-launched attention on integrability in maximally SYM, identified the one- 
loop dilatation operator of scalar gauge-invariant fields (of bare dimension L) with a 
genuine so{6) integrable hamiltonian of a spin chain (with L sites) [2]. In the next 
turmoil, integrable structures were hinted and found in all the sectors of A/" = 4 SYM 
and at all loops (cf. for instance [3]), bearing in mind the convergence of the anomalous 
dimension (weak coupling) expansion. Yet, all these integrable (Bethe, asymptotic) 
scattering equations have the form of a deformation of the one-loop spin chain case, 
plus an additional universal (string) scattering factor; the deformation is such that the 
interaction range increases with the number of loops. Therefore, starting from a certain 
loop order (generically L or higher), they are plagued by the 'wrapping' problem (cf. for 
instance the third of [3]), which was solved only in the SU{2) sector, where a mapping 
to the Hubbard model [1] was possible, however without incorporating the dressing 
factor. In a parallel way, integrability in superstring theory was discovered at classical 
level ItSj and then extended to semiclassical and quantum level. 

In this paper, we want contain ourselves within the (non-compact) sl{2) sector, viz. 
to the local composite operators 



where V is the (symmetrised, traceless) covariant derivative acting in all possible ways 
on the L bosonic fields Z. The spin of these operators is s and L is the so-called 
'twist'. Proper superpositions of operators (11.11) have definite anomalous dimension A 
depending on L, s and the 't Hooft coupling A = Sn^g^: 



where the anomalous part 7 is easily related to the integrable chain energy, E (not to 



(1.1) 



A = L + s + -f{g,s,L) , 



(1.2) 
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be confused with the string energy) via the well-known proportionality 



l{9,s,L) = g^E{g,s,L). 



(1.3) 



Actually, at one loop the integrable problem is equivalent to the analogue in planar one- 
loop QCD for various types of quasi-partonic operators under specific circumstances 



In the context of Bethe Ansatz like equations, a useful tool to perform calculations 
is indeed the so-called non-linear integral equation (NLIE) in its 'excited state' version 
[8]. The NLIE allows to write exact expressions for the eigenvalues of the observables 
for arbitrary values of the system length and of the Bethe root number. Actually, 
this equation turns out to be more efficient for numerical computations as well as for 
analytic evaluations in some particular conditions, e.g. large number of Bethe roots. 
Implementing this idea we previously found and discussed finite size effects to the 
anomalous dimensions in the su{2) sector of A/" = 4 SYM [9l [TOl [11]. In this case, 
the NLIE and the exact expressions for the eigenvalues of the charges have the same 
structure as in models studied in the past. This is a consequence of two simple facts 
concerning this specific case. First one, the scattering matrix between two magnons, 
which appears in the r.h.s. of the Bethe equations, depends only on the difference 
of their rapidities (Bethe roots), provided the so-called dressing factor (cf. below) is 
neglected. Secondly, the Bethe roots completely fill the real axis and only permit the 
presence of a finite number of holes (or complex roots). However, at least one of these 
properties fails when considering other sectors of A/" = 4 SYM, or when the S'-matrix is 
suitably equipped by a string theory CDD factor 0, the dressing term [12] : this is indeed 
the situation in the s/(2) sector. In general, for dealing with this intricated structure of 
the Bethe equations appearing in the whole TV = 4 SYM, we proposed in [13] a path to 
a NLIE substantially different from the original idea of [8]. The different strategy is to 
perform the integrations just on the region (generally intervals of the real axis) in which 
Bethe roots concentrate and, consequently, to avoid the use of the Fourier transform 
in order to write the equation. It follows that this new procedure is effective when the 
magnon scattering matrix has a general dependence on the rapidities and the Bethe 
roots are concentrated on intervals of the real axis or even complex lines, i.e. in all the 
relevant cases of A/" = 4 SYM. As a consequence, we obtain a simplification that we will 
explain here in the sl{2) case: since nonlinear integral terms enjoy an integration just 
on the region where the s Bethe roots actually lie, they become depressed more than 
any inverse power (lns)~", n = 1,2,..., for any fixed value of the twist L. Therefore, 

^Albeit QCD is in the whole not a conformal quantum field theory, it still behaves like one at 
one loop and as far as the anomalous dimensions are concerned. Even under these circumstances, 
integrability is not complete since it requires additional constraints as, for instance, aligned helicities 
of the partonic degrees of freedom. 

^Although originally thought of as a correction to the scattering coming from string effects, in the 
end it revealed its efi'ects already at four loops. 
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stated in advance the logarithmic scahng 

j{g,s,L) = f{g)lns + 0{s'), (1.4) 

after the very important BES' paper [12] on f{g) (cf. also e.g. [H] for preceding 
literature), we are left with a Linear Integral Equation (LIE) which allows us to compute 
the sub-leading corrections, 0{s^), of the conformal dimensions and, of the other 
charges, both the leading and the sub leading terms. In this respect, our approach is 
clearly different from that of [16], which uses the full real axis NLIE presentation by 
[8], because this needs to take into account and to evaluate the non-linear integrals (on 
the whole real axis) as well. 

Very interestingly, our LIE does not differ from the BES one [12] , but for the inho- 
mogeneous part, which consists in an integral on the one loop root density and a hole 
depending term (apart from a known function) 0. As the equation is non-perturbative 
(i.e. for any g) and linear and it drives the cross-over between weak (small g) to strong 
(large g) coupling in a rather intelligible way. Nevertheless, for evaluating and check- 
ing dominant string effects (like for instance the dressing phase) the aforementioned 
wrapping effects ought to be negligible or known. 

Eventually, these facts have furnished us the stimulus to investigate the next-to- 
leading-order (nlo) term - although coming from an asymptotic Bethe Ansatz - in that 
the leading order f{g) has been conjectured to be independent of L or universal [12] , 
after the one loop proof by 

In this paper we will give an explicit application of this new type of (N)LIE to the 
(L, s)-vacua of the s/(2) sector with the following plan. In Section 2, we will outline 
the formalism suitable for writing the NLIE on an interval in the most general case. In 
Section 3 we will apply this technique to the case of the spin —1/2 XXX chain, which 
describes the one loop s/(2) sector of A/" = 4 SYM and in Section 4 we will discuss 
the many loop case and end up with our main object, the linear integral equation. By 
means of the latter, we will compute up to three loops in the 't Hooft coupling the 
leading and sub-leading corrections of the eigenvalues of the conserved charges in the 
large s limit, and check our results vs. the anomalous dimensions in [16] (for what 
concerns the other charges, our findings are new at the best of our knowledge). 

^Then, we could also consider the limit L = jlns oo for fixed j, which is indeed the rele- 
vant scaling of this theory [15j and entails an improvement of the previous formula into "/{g, s, L) = 
f{g,j)\ns + 0{{\n s)~°°) [16]; but this would be the subject of some future publications. 

■*This structure encourages us to proceed further in the direction of the fixed j expansion and 
suggests the preservation of itself: but this way may only be the topic of future publications. 
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2 A new approach: the NLIE on the interval (s) 



In almost all the cases considered up to now, the NLIE was written for counting func- 
tions defined as 

s 

Z{u) = <t>{u)-Y,<P{u-Uk), (2.1) 

k=l 

and when the Bethe roots distribute on the real axis, allowing the presence of only a 
finite number of holes and possibly complex roots. Even if this case is relevant for the 
study of the fundamental state and the first excitations of many models, it does not 
cover many of the Bethe Ansatz systems proposed in the context of A/" = 4 SYM. 

For this reason we want to write the NLIE (and the expression for the eigenvalues of 
the observables in terms of its solution) for the more general case in which the counting 
function is defined as 

s 

Z(m) = $(«)- ^0(M,Mfc), (2.2) 

k=l 

(i.e. the function 0(x, y) does not depend only on the difference x — y: this happens, 
for instance, when the dressing factor is present). We suppose also that the s Bethe 
roots {uk}k=i,...,s are concentrated in an interval [A, B] of the real axis and that a 
finite number of holes is present. We call the holes (in number Hi) lying inside the 
interval and the holes (in number Ho) lying outside the interval. This particular 
distribution of roots is peculiar, for instance, of states in the sl{2) sector of A/" = 4 
SYM. 

On this state we consider a sum over the Bethe roots {uk}k=i,...,s of a function 
(observable) 0{u) analytic in a strip around the real axis. Without loss of generality, 
we can put ourselves in the case in which on both roots and holes the condition e*^'-"'''' = 
giZ{n^ ) = —1 holds. In this case, this sum can be written [8] as 

2mJ20iu,)=\nn / rf.O(. - .e)^--^A__J+ (2.3) 

k=l '-"'A "T 



+ / duO{u + ie) — ^— — ^ -27rz> O 

''^ -I h=l 



Supposing Z' (u) < and supposing that the values Z{A) and Z{B) are known, we can 



^The case in which the Bethe roots are concentrated on a finite number of intervals on the real axis 
follows straightforwardly from the results of this Section. Moreover, even the case when the roots lie 
on complex lines can be treated as follows. 
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rearrange this expression as follows, 



= - f ^0{v)Z'{v)+ [ ^0(^)imln[l + e'^(-°)]- 
^ Ja Ja ^ ""^ 

- f^0i4^) = -^[0iB)ZiB)-0iA)ZiA)]+ (2.4) 



+ ^ {0(5)Imln [1 + e'^(^)] - 0(A)Imln [l + e^^^^)] } + 



+ 



27r 



0'(t;)Z(t;) - 2 / ^0'(t;)lmln [l + e^^(^-^°)] - f]0(M«) . 

h=i 



In brief, what we are doing is to evaluate a sum on the Bethe roots by integrating just on 
the interval containing them. Therefore, this method is alternative and complementary 
to the idea proposed in the first of [8] which consists in first integrating on all the real 
axis and then subtracting the contributions coming from the real holes. The convenience 
of this new method is that the non-linear terms present in (12.41) are strongly suppressed 
in the limit of large number of Bethe roots. We will come back on this point in the 
next sections, when we shall apply this technique to the s/(2) sector of A/" = 4 SYM. 

We now use ( 12. 4p in the sum over the Bethe roots appearing in the definition (12. 2p 
and obtain the following equation 



(i) 

h=l 



Z{u) = $(«)+ [ ^<P{n,v)^Z{v) + f^ 
■^^ ^ ^ h=l 

_ 2 r — 0(n,t;)— Im In [1 + 6*^(^-^0)1 = 
Ja 2n dv ^ 

+ 2 Tjl^-^^ii.i'jlmlnh+e'^'"-"'] , (2.5) 



where 



f{u) = ^{u) + ^[ct>{u,B)Z{B)-ct>{u,A)Z{A)] + Y.4>{u,uf)- 
Ztt ^ — ' 

h=l 

- - {(f){u, 5)Imln [1 + e*^(^)] - (j){u, A)Imln [l + e*^^^)] } . (2.6) 



We can now write a NLIE for the counting function by inserting in an iterative way 
( 12. 5p for Z in the right hand side of the same equation. Using the notation 

i^*f)iu)= r dv^iu,v)fiv), (2.7) 
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eventually we gain the NLIE in the final form 

Z{u) = F{u) + 2{G^L){u), {2.t 

where 



(2.9) 



[U,V) 

k=l k=2 



We used the simplified notations 

L(u) = Imlnh +e'^("-'°)l , ip{u,v) = ^^(f){u,v) . (2.10) 

More explicitly, 

F{u) = /(u) + V(-l)M dviip{u,vi) dv2v{vi,V2) ... 

k=i 

dvk(p{vk-i,Vk)f{vk), (2.11) 

'A 

pB pB 

G{u,v) = v9(M,t;) + y'(-l)^ / dvoip{u,Vo) dvii^{vo,vi) . . . 

k=i -^^ -^^ 

-B 

dvk-i 'fivk-2, Vk-i) (fivk-i, v) . (2.12) 

These expressions are quite formal and very difficult to handle. More easily, both the 
forcing term F{u) and the kernel G(u, v) may be found by solving respectively the linear 
integral equation, which easily follow from ( 12.111 12.12p . 

F{u) = f{u)- f dvip{u,v)F{v), (2.13) 

J A 

.B 

G{u,v) = ip{u,v) — / dwip{u,w)G{w,v) , (2-14) 



(2.15) 



and are linked the one to the other via 



F(m) = /(«)- / dvGiu,v)fiv). (2.16) 

J A 

Eventually, inserting (12. 8p in (12.41) we obtain an expression for the eigenvalues of an 
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observable as 

= -^[0{B)Z{B)-0{A)Z{A)] + 



k=l 



271 

+ - {0(5)Imln [1 + e^^(^)] - 0(A)Imln [l + e'^^^^] ] + 



+ / ^0\v)F{v)-Y,0{uf)+ (2.17) 

pB J rB 

+ 2/ —0\v) dw[G{v,w)-5{v-w)]\m\n[l + e'^^'"-'''^] . 

We remark that all the already known NLIEs can be reproduced in this way, without 
Fourier transforming. In this sense the method sketched in this section is more general. 
It seems natural to use formulae (12. 8[ 12. 9p and (12.171) in order to write, respectively, the 
NLIE and the eigenvalues of the observables on states appearing in models relevant for 
A/" = 4 SYM. In the next two sections we will apply this new method to the widely 
studied |3l [13 [I2] sl{2) sector of the theory. 



3 An example: the XXX_i spin chain 

We want to apply the techniques developed in the previous Section to the XXX_i 
spin chain, as the latter clearly gives a representation of the s/(2) sector at one loop, 
as stated in the 'Introduction'. 

The spectrum of this non-compact spin chain may be described by the Bethe equa- 
tions, 

where we have indicated with L the length of the chain and s the number of Bethe roots. 
In this case the s Bethe roots concentrate in an the interval of the real axis symmetric 
with respect to zero. L holes are present: two holes uf^\ h = 1,2, lie outside. Hi = L — 2 
holes u^i^\ h = 1, . . . , L — 2, lie inside this interval. Let us define the counting function 
(for reasons that will be clear in the following, we will put an index to all the functions, 
e.g. Z, F, G, related to the XXX_i spin chain) as 

.„,.)..,m^)-.p„(i±^). ,3.) 
From the Bethe equations (13.11) . we have the conditions 



ainr^^^]-ty\nr^^^^]=27rn,, E Z . (3.3) 
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^ln(^^^-^lnf^^^)=7r, ^>0, (3.4) 



Using the property 



X + J V^^ + ^- 

we have that, on both Bethe roots Uk and holes Uh, 

Zo{uk) = n{2nk - L - s + 1) , Zo{uh) = n{2nh - L - s + 1) . (3.5) 

We choose L + s even, in such a way that e*^''"*^ = e*^^"h'^ = e*^^"h'^ = —1, as 
in the previous Sectioiij. If we now go from the smallest to the biggest root of the 
interval (passing also L — 2 internal holes) the (decreasing) counting function Zq varies 
of —271(5 + L — 3). Since Zq is an odd function, this means that 

Zo{us) = -rt{s + L - 3) , Zoiui"^) = -7r(s + L - 1) , (3.6) 

where Ug is the biggest (most positive) root of the interval and u^°^ is the positive hole 
outside the interval. We can choose the separator bo, which defines the interval of 
integration in our formulae, i.e. A = —bo , B = bo, such that 

Zoi-bo) = -Zoibo) = 7ris + L-2). (3.7) 

With this position, the relevant functions defined in the previous section take the form 

^oi^) = — 2L arctan2'U , </)o('U, f ) = 2 arctan('U — t>) , (fo{u,v) = —, 

n 1 + [u — vY 

fo{u) = — 2Larctan2M — (s + — 2)[arctan(n — 60) + ai'ctan(n + 60)] + (3.8) 

L-2 

+ 2 arctan('U — ■uj^*'') . 

h=l 

Since s + L is even, we have that Imln [l + e*^o(±feo)j _ Jj^ addition, one easily shows 
that in this case the functions Fo and Go enjoy the parity properties 

Fo(m) = -Fo(-m) , Goiu,v) = -Goi-u,v)=Goi-u,~v). (3.9) 

They may be determined by solving the linear equations 

Fq{u) = —2L arctan2M — {s + L — 2)[arctan('U — 60) + arctan('U + 60)] + 

+ 2^arctan(n-n«)+ / --— r^Fo{v) , (3.10) 



h=l 



1 1 r^" dw 



Go{u,v) = — — — — — — -^Go{w,v). (3.11) 

TT 1 + {U - V)^ J_hg IT 1 + [U- WY 

At the leading s — 00 order these equations become simpler and therefore we will start 
with this case. 



^The number of derivatives s has to be even and in any case, there is no loss of generahty, in that, 
if L + s is odd, the only modification in all the formulae is the replacement of the logarithmic indicator 



Lq{u) by Imln [l 
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3.1 Determination of Fq when s +00 

Let us consider equation (I3.10p for the function Fo{u). It is convenient to rescale the 
variable u = us and the extreme of the interval 60 = ^os, in view of the limit s — > +00: 

Fo{us) = —2L arctan 2us — (s + L — 2) [arctan('us — 60s) + arctan('us + bos)] + 

+ 2 Varctan(us-4'^) + / dv-- — -Fo(t;s) . (3.12) 

Using the asymptotic expansions (0(s~") means "terms of order s~"") 

71 1 

arctan-us = — sgnu — ^ — \- 0(s^^) (3.13) 
2 us 

-T—^^—V2 = S{u-v)--^-^ + 0{s-'), (3.14) 
TT 1 + — v)"^ Sir duu — V 



msi 



de the integral and the fact that 



lim = 0, (3.15) 



we obtain in this approximation 



^ , , , , 4:- L s + L-2 2u 

Fo{us) = -27rsgn(n) + ^— + — — ^ + Fo{us) + 

us S — On 











dv — j 








\dv u — V J 



+ j dv— [ —P- — - ] Fo{vs) . (3.16) 

Integrating by part gives 

A-L s + L-2 2u Foi-bos) f 1 , 1 
= — 27rsgn(u) H ^ h 



us s u^ — b^ SIT \u + bo u — bo 

dv—iP—]—Foivs). (3.17) 

Sir \ u — V J dv 

Now, condition (13.71) at large s imposes on the function Fo the constraint 

Foihs) = -Foi-bos) = -7i{s + L - 2) + o(s°) , (3.18) 

where o(s°) means "terms of order smaller than . Therefore, equation fl3.17p simpli- 
fies into 

A — T r^o 1/ T \ d 

= -27rsgn(u) + / dv— ( P ] Fo{vs) . (3.19) 

us J__lg STl \ u — V J dv 



'^This is indeed a subtle exchange of two limits. 



10 



This equation can be solved by finite Hilbert transform tecliniques. Its solution reads 

2 



1 d 
27CS du 



Foius) 



In 



&0 + 



TT 



U 



2 - ^ 



(3.20) 



Foius) 



u 



Ah^s arcsin 2us In 



^0 



TT ( 2 - ^ ) sgnw . 



Remember, however, that also 6o depends on s through (13. 7p . which at large s implies 
fl3.18p - a normalization condition for f l3.20p in this approximation 



6o 



du 



bo 



-In 

TT 



&0 + 



2 - ^ 

^ 2 



1 + 



L-2 



(3.21) 



From this equation we deduce, after integration. 



2bo 



2-# 



L-2 



J ^ ( L 



(3.22) 



Inserting this expansion in (13.20p . we obtain the following behaviour of Fq and its 
derivative when s ^ oo: 



1 d 



2tis du 



Fn(us) 



1 



In 



TT 



Foius) 



2 ' V 4 ""^ 



U 



arcsin 2^ + m In 



+ 



1 



2 - 



TTS /l_-2 S 



-2s 



TT I 2 - — 1 sgnti . 



u 



2 arcsin 2u + 



(3.23) 



3.2 The leading order density equation from the NLIE 

Once we have determined the functions -Fo('^^) ^^^d G'o(m, f), as well as the value of the 
extreme ho-, one can write, according to the general formula 02.80 the nonlinear integral 
equation on a finite interval satisfied by the counting function of the XXX_ii2 spin 
chain. It is not difficult to relate such NLIE with the linear equation [17] satisfied by 
the density of roots in the limit s —>■ +oo. In this respect, it is convenient to use the 
first equation in the chain (12. 5p 



Zn(u) = — 2Larctan2M+ / — 2 arctan(M — t;) — [Zo(f ) — 2Lo(f )] + 



-bo 



27r 



dv 



L-2 



+ 2 arctan(M — -uj^*^) 



(3.24) 



h=l 
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where Lq(v) = Imln[l + e«^o{^'-«o)j^ j^^ before, we rescale the variable u = us and 
the extreme of the interval Bq = BqS and then we let s +00. Using the asymptotic 
expansion 

TT 1 

arctan-us = — sgn-u — ^ — h 0(s~^) (3.25) 
2 us 

inside the integral and the expansion fl3.22p for 60, we obtain, in this approximation, 
the equation 



Zo{us) 



' 2 dv 

—TcLsgnu + I — 



1 27r 

2 



vrsgn M — VI — 



siu — V) 



-^[Zoivs)~2Loivs)] + 
dv 



L-2 



(3.26) 



h=l 



Performing the integral involving the sgn function and using the fact that, when s 



-00, u 



0, we are left with 



2Ln(us) 



f 2 dv 


2 


/ 1 2^ 

2 


s{u — v) 



^[Zoivs)-2Lo{vs)]. (3.27) 



As the term involving Lq{vs) approaches zero (cf also the Appendix), we can neglect it 
and obtain 



Defining the density 



= — 27rsgnM 



2 dv 



d 



27r s{u — v) dv 



Znivs) . 



-Z;o 



1 d 



-Zo{us) 



(3.28) 



(3.29) 



27TS du 

which, because of the condition ^o(^o) ~ Zq^—Bq) = —2n{s + L — 2), is normahsed as 



dupo{u) = 1 + o(s°) 



(3.30) 



our equation is written as, 

= -27rsgnM + 2 



1 

2 _ 1 _ _ 
dv- -po{v) , 

1 u — V 

2 



(3.31) 



which coincides with (52) of [17\. Therefore, the NLIE on a finite interval for the 
XXX^i/2 spin chain links itself in a simple way to the linear equation for the density 
of roots in the limit s —>■ +00. 
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3.3 Determination of Go when s ^ oo 

On the other hand, Go{u,v) satisfies the integral equation 



Go(u,v) 



1 



1 



TT 1 + (m — v)"^ 

which, in terms of rescaled variables reads as 



dw 1 
vr 1 + (m — wY 



Go{w,v) 



(3.32) 



-bo 



Gn(uS, VS) 



1 



+ 



dw s 



Ti 1 + s'^{u — vY 

Using the expansions f l3.13p . we obtain the equation 

1 \ f d 1 \ 

= --5{u -v) + — dw( —P Go 

s ns J \dw u — w J 



Go{sw, sv) 



(3.33) 



sw, sv) + o{s~ 



We now define the order s° antisymmetric combination 

Xo{u, v) = Go{su, sv) — Go{su, —sv) + o(s°) , 
which satisfies the equation 



(3.34) 



(3.35) 



1 1 



6{u + v) -6{u-v) + -- ^Xo -,v 



TT U 



1 1 



vr u 



'2' 



V - 



1 



dw P 



2 



1 



u — w J dw 



d 



-Xn(w,v) = 0. 



(3.36) 



We have set bo = 1/2, since Xq takes into account only the leading 0(s°) contribution. 
The solution to this equation is given by 



11 11 



^Xo{u,v) 



+ 2 

Integrating this function we may write down 
1 , V + u 



TT V — U TX V + U 



6 { \ + 6 { u + - 



\-v^ 



Xo(u,v) 



In 



TT 



V — U 



-7 — \ / i — u"^ — vu + 



i-v^Jl-u^ + vu + \ 



+ 



X,{-,v 



2v 



(3.37) 



arcsin 2u + 



+ 



sgn I w — - I + sgn { u + - 



Xui-,v 



(3.38) 



where we took into account the fact that Xq{u,v) = —Xq{—u,v) in order to fix the 
undetermined function of v. 
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We also define the order s° symmetric combination 

Yq{u, v) = Go{su, sv) + Go{su, —sv) + o(s°) , (3.39) 
which satisfies the equation 

-i9(u-v) +ie(-u-v) = — dw ( P-^]yo(w,v) . (3.40) 

«vr J_i \ w — u J 



This equation can be solved as 



YJu.v) = — [' dw(p-^^\-t9(w~v)+i9(-w~v)]J4 — ^= (3.41) 




We remark that this solution satisfies the parity properties Yo{u,v) = Yo{u,—v) = 
Y{—u,v), as it follows also from properties (13.91) . As boundary condition, we find 

Yo (±1, v] = Yo (u, ±^1=0. (3.42) 



. 2' y V 2^ 
3.4 Evaluation of the charges when s ^ oo 

Using (12.1 7p we want to compute the eigenvalues of the energy and of all the charges 
when s oo. We are interested in the leading terms in the large s expansion, i.e. in 
the terms proportional to Ins and to s°. We have to specialise formula (I2.17P to the 
case in which 

A = -bo, B = bo, Zo{-bo) = -Zo{bo)=7T{s + L-2), (3.43) 



0{v) = qr{v) 



I 



r — 1 



1 1 



(3.44) 



For parity reasons, the eigenvalues of the charges Qr-, with r odd, are zero. Therefore, we 
restrict to even r. One easily sees also that the first two lines of (12.171) give contributions 
at most 0{s^^). More importantly, we have strong evidence from numerical simulations 
that the nonlinear terms in the fourth line of (12.171) go to zero as s ^ oo. This peculiar 
behaviour of the nonlinear terms is due to the fact (see for example the discussion in 
Appendix A of [12]) that in our approach we are integrating only on the interval in 
which the Bethe roots and the internal holes are present and not - as in the approach 
of [8] used by [16] - on the whole real line. It follows that, differently from at 
least up to the order o{s^) at which non-linear terms in the fourth line of (I2.17P start 
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contributing, the eigenvalues of the charges are given by only the linear term in the 
third line (cf also Appendix) 



Qr 



h=l 



bo 



2n dv 



r — 1 



(3.45) 



which involves Fq, i.e. the solution of the linear integral equation (13.121) . 

One could now insert for Fq the solution at large s given by formula (13.201) . However, 
we will now show that this procedure is accurate only for the energy (r = 2) and 
only for the (coefficient of the) leading Ins term. Indeed, let us use (I3.20p . In 

fact, we obtain in this approximation 

dv 2i \ 1 1 ] [. fbo + ^/bl-v^^ ^ 



Qr 



-bo 



27rr - 1 



L 



r-l 



In 



- n{2--]6{v 



5 Abot 
dx 



-^(L-2)2^ = 
1 2-(-l)i-i L 

— ^/^o(-) - -j^-Yr2 



r - 1) {2boX + i) 



identifying the Korchemsky density of Bethe roots (s — > oo) [18] as 

Po{x) = - In 



X 



(3.46) 



(3.47) 



Introducing the resolvent 

G{x) 

we see that 



dy 



Po{y) 



i In 



VI - 4x2 + 1 



1 y-x 



Qr 



2(26o)^~" d"--"^ 



;G{x)l 



I (r-l)!dx-2"-^l-=V4''o 
Explicit computation would give, for r = 2, 



VI - 4x2 - 1 

2^(-l)i-i L 
" (r-l) 2 



+ o(s°) . 



Q^ = E = 41ns + 41n2 -2L + o(s°) , 



and, for r > 4, 

2^(-l)i- 



Q 



_ 2'{-iy-iL 
(r - l)(r - 2) (r-l) 2 



o[s 



2^(-l)i- 



(r - l)(r - 2) 



(3.48) 



(3.49) 



(3.50) 



(3.51) 



15 



One can check for L = 2, 3 that the order terms in the energy and in the higher 
charges are not the correct ones (cf. for instance [HI [T71 [IHl [20]). Therefore, the 
large spin solution for Fq , fl3.20p , is a good approximation only if we are interested in the 
leading O(lns) term of the energy, although Fq would be exact up to the 0(s°) order, 
i.e. neglecting (the going to zero) o(s°) terms. Yet, this failure reflects the subtlety of 
considering the s — > +oo limit of the integral equation (13.101) . 

Nevertheless, a sufficiently accurate approximation for Fq, efficient when neglecting 
in the charges the o(s°) terms |§, comes out by solving the one loop density equation, 
i.e. the derivative of (both members of) equation (I3.10p . upon approximating m^*'' = 0, 
by means of Fourier transform technique 



e 2 cos, ,1.^/ V ^ ; o - r, 

ikFQ{k) = -4n^ ' ^ +2n{L-2) ^-(47rln2)(5(fc)+o(s°) . (3.52) 

2 sinh 2 sinh ^ 

Using (13.521) and the Fourier transform of the charges functions Qrik), the eigenvalues 
of Qr at order In s and are given by 

—q,(k)tkFo{k) - ^-y^iL - 2)r + o(s°) . (3.53) 

From (13.531) we obtain the correct results for the eigenvalues of the charges 
Q2 = E = 41ns + 47£;-4(L-2)ln2 + o(s°), 

Qr = ^tiL!^l!:illl[(2-2'-i)L-2(l-2'-i)] + o(s°), r>4, (3.54) 
r — 1 

where 'Je is the Euler-Mascheroni constant and ({x) is the Riemann zeta function. 



4 All loops 

Let us now consider the Bethe ansatz like equations 
where we used the notations 



A = Svr^/ , (4.2) 

A being the 't Hooft coupling. It is believed that configurations of Bethe roots, i.e. so- 
lutions of (14.11) . and the corresponding eigenvalues of the energy are related respectively 

^More precisely, the first order we are neglecting is 0(1/ Ins) and comes from the approximation 

. . (i) 

of all the internal hole positions uj^ =0, which is strictly valid only if s = +00. 



x^ = x^{uk) = x{uk±i/2) , x{u) = - 



1 + A 1-^ 



2g^ 
2 



u 
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to composite operators and their anomalous dimensions in the si (2) sector of jV' = 4 
SYM. This correspondence, however, breaks at the so-called wrapping order, i.e. at 
order g"^^'"^, and higher. Therefore, results in this section are relevant for jV' = 4 SYM 
only until the order g^^^"^. 

In a fashion similar to the one loop case, Bethe roots concentrate in an interval 
[—b, b] of the real axis. Inside this interval, L — 2 holes arc present, while outside it two 
external holes lie. We use for them the same notations as in the one loop case. 

The counting function is 



Z{u) 



—2L arctan 2u — iL In 




) 



2 arctan(M — Uj) + 



s 



+ 




(4.3) 



where the so-called dressing factor is given by 



oo 



oo 



■+l+2u{g)[qr{Uk)qr+l+2u{Uj) - qr{Uj)qr+i+2u{Uk)] , 



(4.4) 



r=2 i/=0 



the functions Pr,r+i+2u{g) being 





(4.5) 



and qr{u) being the density of the r-th charge 




(4.6) 
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H 



4.1 Equations for F, Fq and F 

The counting function (14 .Sp can be treated according to the general formahsm given in 
Section 2, provided we choose A = —b, B = b and 



1 + 



$(u) 
(j){u,v) 



—2L arctan 2u — iL In 



2x-(uY 



1 + 



2 arctan(M — v) 



2i 



In 



1 3L 



2x+(u)x-{v) 



1 - 



+ i6[u, v) 



2x (n)a;+(i;) 
1 + ^ 



—2L arctan 2u — iL In 



2x-{uy 



1 + 



+ 



2x+{uy 



L-2 



/l=l 



arctan(-u — u^^'' (g)) — i In 

arctan(M — b) — i\n 
1 - 



1 



+ -Z{b) 



TT 



+ 



2x+(«)x-(n«(g)) 



1 - 



2x-Hx+(n«(g)) 



2x+{u)x-(b) 



1 - 



+ e(?/,6) 



2x-(u)2:+(6) 



+ 



arctan(u + 6) — z In 



1 - 



2x+{u)x-{-b) 
2x-{u)x+{-b) 



+ 9{u, -b) 



(4.7) 



with the exphcit (^-dependence of the internal all loop holes u^^{g). We have sup- 
posed that Imln[l + e'^^^^^] = 0, as in the one loop case. It follows that the function 
F{u) entering the NLIE satisfies the linear equation (12.131) with the function f{u) and 
ip{u,v) = 2^^0('W)'^^) obtained from (14.71) . We can split F{u) into its one loop contri- 
bution Fq and its higher loop contribution F^{u): 



F{u) = Fo{u) + F"{u) . 
Of course, the one loop contribution satisfies the LIE (13.101) 

Fq{u) = —2L arctan 2u — (s -|- L — 2) [arctan(u — 6o) + arctan(u -|- 6o)] + 

L~2 „fc,) ^ ^ 

+ 2Varctan(M-M«(0))+ / dv-—- r^F^iv) , 



(4.9) 



where ^^^^''(O) are indeed the internal one loop holes (as g = value of the (internal) 
all loop holes). On the contrary, the LIE obeyed by the higher loop F^{u) contains 
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additionally the (^-depending holes in the form 



1 



U 



—iL In 



+ i arctan('U — u^^\g)) 



2x-{uY 



2x+{uY 

i arctanf-u 



L-2 / 1 



2x+(u)x-(u^^\g)) 

2X-(«)x+(«W(g)) 



In 



1 - 



2x^(u)x-{h) 



In 



<(0)) 
1 



2x+{n)x-{-fe) 



2x-(«)x+(fe) 



1 



+ i9{u,b) + i9{u,-b) 



+ 



2x-{u)x+{-b) 

+ — Z(6)[arctan('u — b) + arctan(-u + b)] Z(6o)[arctan(M — bo) + arctan('u + bo)] + 

TT ~ 



TT 



+ 



+ 



dv 



F^'iv) 



-b 



TT 1 



dv 



-bo 



dv 



TT 1 



v)' 



:Fn(v) 



bo 



^ dv 1 

TT 1 + (m — t>)2 



(if 



In 



2x'^{u)x {v) 



d 



+ i—6{u,v) 
dv 



[Fo{v) + F^{v)] 



Fo{v) + 
(4.10) 



2x~ {u)x+ (v) 

Now, we proceed to the limit s +oo keeping, of course, only the non-vanishing 
terms. Upon treating by parts the last integral, the term in the square brackets in the 
third line of (I4.1UI) gets multiplied by the factor Z{b) — F(b) = o(s°)§|. Therefore it can 



be neglected. Furthermore, since uf'{g) = 0(5*^), we can set everywhere u^^id) ~ ^^ 
Moreover, as s ^ +00, one has that Z(b) = 0{s), 2'o(feo) = 0{s), but also that 
arctan(M — b) + arctan(M + b) = 0(s^^), arctan(M — b^) + arctan(M + bo) = 0(s^^). So, 
the terms in the fourth line can be neglected as well. Since b — bo is 0(s°) and -Fo(^) is 
0{s), the second and the third integral in the fifth line can also be removed. For the 
same reason, in the integrals in the sixth line involving Fo we can set the extremes to 
their one loop value =f6o- Finally, in the integrals involving F^{v) we may replace the 
extreme, b, - which in principle ia an unknown of the integral equation fl4.10p as well 
- with +00, since as f — 00 F^{v) becomes constant, while the rest of the integrands 
vanishes as or faster. Rearranging the terms in the right hand side, we have 



F^iu) 



1 + 



-iLln 



2x-{uy 



2i{L - 2) 



I 

TT 



pbo 




/ dv 


In 1^ 


'-bo 





2x+{uy 
1 - 



In 



1 - 



+ 



dv 



vr 1 



dv 



2x+{u)x-(v) 
2x 



+ i9{u, v) 



2x+(u)x-(0) 
2x-(u)x+(0) 

^Foiv) + 
dv 



+ ie{u,0) 



[V - 



[4.11] 



In 



1 



2x+(u)x~ (v) 
9^ 

2x (u)x''^{v) 



+ i9{u, v) 



dv 



^This is certainly true at one loop, as follows by comparing (|3.7[) and (|3.18p . We suppose that it 
stays true also when considering all the loop corrections. 
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With such rearrangement, we have collected in the first two lines of the right hand side 
the forcing terms of our equation. In the last two lines, we have the integral terms, i.e. 
terms involving integral of a kernel function with the unknown . 

Several comments are now in order. The LIE (14. lip constraining the density of roots 
— may be thought of as an improvement of the BES equation [12], which 

in fact takes into account not only the leading O(lns) term, but also the sub leading 
0{s^) (constant) corrections. This equation is exact and linear and therefore drives the 
interpolation between weak (small g) to strong (large g) coupling in a non-perturbative 
way. In particular, it may be useful for studying the strong coupling regime where 
the string effects, and in particular the dressing phase, dominate. However, in view of 
possible comparisons and checks vs. string results, one should fucus the attention in 
cases where the aforementioned wrapping effects are negligible or known. 

As in the one loop case, in order to compute the energy and the eigenvalues of 
the charges in the limit s ^ oo and including constant terms, in our framework it is 
sufficient to consider the third line of the general formula (I2.17p . Explicitly, 



E{g,s) 



^» dv d 
-bo 271" dv 
dv d 
, 271 dv 



x+{v) x~{v) 
i i 



x^iv] X [v] 



Fo{v) - (L - 2) 



x+{0) x-(0) 



+ 



F^'iv) + o(s^) . 



(4.12) 



This means that, in contrast to the approach of [TB], where the [S] method is used, we 
have to cope only with linear equations. In respect to (14.120 and for future convenience, 
we introduce the function 



h{g,s) 



bo 



-bo 



dv 
2^ 



^X^[V) X [V 

dv 2i d ^ , , ,^ , 1 
Fo(w) + 2i(L - 2) 



27r X iv) dv 
in such a way that the energy reads 

, 2^ 



x-(0) ' 



x+(0) x-(0) 



E{g,s) = h{g,s) 



x+(f) x^(v) 



-^F^(t;) + o(s°). 
dv 



(4.13) 



(4.14) 



Direct calculations show that 



h{g,s) = 41ns + 47£;-4(L-2)ln2 + /i(c/), (4.15) 
hig) = 2/[3LC(3)-7C(3)] + (7^(62 - 30L)C(5) + O(/) + o(s0). 

In addition, we remark that equation (14. lip has the same kernel as the BES equation. 
They differ only in their forcing terms, since in the BES case, in contrast to (14. lip , the 
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forcing terms are simply 



2 r 



271 

2^(2 Ins 



1 1 

+ 



1 



2i 



TX 



u=0 



bo 



dv- 



d 



-bo 



V — ^dv 



+ 



+ 8 In S ^ f32,3+2.{g)q3+2u{u) + 0(s°) 



(4.16) 



u=0 



where we used the one loop results at order Ins: 



bo 



bo 



dv 



d 



27rv + ^ dv 



Foiv) 



bo 



bo 



dv 1 d 
271 V — ^ dv 



Fo(t;) = 2ilns + 0(s°). 



(4.17) 



Our forcing terms differ from the ones in the BES equation for two reasons. Firstly, in 
addition to the BES terms, we have genuine new terms coming from the first line of 
the right hand side of (14. lip . If we expand them in powers of g^, we see that they are 
structurally different from the BES forcing terms, with the exception of the term 



2i{L-2) 



g 



2x+{u)x-{Q) 2x-{u)x+{Q)_ 



2(L- 2)^/32,3+2.(^7)93 



u=0 



^g_ 

2tt 



2 r 



1 1 

+ 



x^[u) 



x~{u) 



2i 



TT 



^p2,3+2u{g)q3+2p 



[U 



u=0 



27i{L - 2) 
x-{0) 



(4.18) 



which is proportional to the BES forcing terms. Secondly, for what concerns the terms 
in the second line of the right hand side of (14.111) , in their expansions in powers of g"^ all 
the terms have to be kept, since we want to be precise in s up to the order s°. However, 
some of these terms, namely 



271 



1 



1 



X+('U) X {u) 



2i 



vr 



^(i2,3+2u{g)q3+2v{u) 



v=0 



dv- 



d 



-bo X i'")dv 



Fo{v) 



as functions of u, are also proportional to the BES ones. We conclude that in our 
equation (14.111) a part of the forcing terms. 



2n 

bo 



x^(u) 



X (m) 



97 

u=0 



2,3+2i^ 



{g)q3+2u{u) 



271 



1 1 

+ — r 



x^iu] 



d ^ , , 27r(L-2) 
X (0) 

2^^^ 

H / ^ P2,3+2iy 

u=0 



X [U] 



vr 



{g)q3+2u{u)}-h{g,s), (4.19) 



is proportional to the forcing terms in the BES equation (the forcing term in the BES 
equation is given by (14.191) in which h{g,s) is replaced by 41ns). Therefore we are 
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allowed to say that the solution of our equation is the sum of a function proportional 
to the solution of the BES equation and an unknown function, 

h{g,s) 



-oo, F"{u 



41ns 



(4.20) 



This means that in the expression for the energy E{g,s) we expect that 
E{g,s) = h{g,s) + h{g,s)[f{g)-A]+E^^'^'^{g,s) + o{s') 

= l%,,)/(^) + E-*™(^,s) + o(s°), 



(4.21) 



where f{g) is the universal scaling function of A/" = 4 SYM and E'^^^^"'{g, s) indicates 
contributions coming from F^^^™{u). 

We can find the structure fl4.2ip for the energy E{g, s), by performing a hrute force 
perturbative expansion in order to solve equation ( 14. lip and to compute the energy. 
Let us define 

F^{u) = g'F^iu) + g'F^iu) + 0(/) . (4.22) 
At the order g"^ we have the equation 



^ ^ ' 2i 



1 



+ 



dv 



1 



- 2(L-2 



1 



1 1 \ 1 1 /"^ , 1 ^ , , 

dv T^F^{v) 



u + 



u 



2m n — f V + ^dv 



^ 1 d 
dv -—Fn{v) + o(s^) 



2in u + ^ J^b V — ^dv 



(4.23) 



2^-0 - 2 

We can now use the one loop results and then pass to the Fourier transform 



![:^e-^+e-l'=lFf(A;) 



IM 

2 



- 47r2[lns + 7£; - (L - 2) ln2]sgn(A;)e 
Solving this equation and going back to the coordinate space we reach 



(4.24) 



{u) = 27r[ln s + 7^ - (L - 2) In 2] tanh un + — 



In a similar fashion, one computes F^{u). We omit the details and give only the final 
result: 



Ail 

16 du^ 



12^ 'du 



+ 



+ -[lns + 7i5-(L-2)ln2] — tanhvTM 



TT |y [Ins + 7£; - (L - 2) ln2] + 7C(3) - 2LC(3)| tanhyrti . (4.25) 
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Therefore, we are allowed to write that 



7T . (P 



F^{u) = [\ns + 'jE — {L — 2) In 2] ( tanh nu H — 9^'i~2 ^^^^ ™ 9^ tanh ttu 

2 du 3 



7r(7C(3) - 2LC(3))/tanh7rM + 



+ 



d 

du 



1 



Tpy- — iu]+'il)\- + iu 



1 



+ 



Ail 

16 c/m^ 



1 



(4.26) 



After obtaining these results we can evaluate the energy E{g, s) up to the order and 
at the orders Ins and s". From the general formula (12. 171) we may write that 



E{g,s) 



dv d 
2tt dv 



[F,{v)+F"{v)]-{L-2) 



X+(f) x~{v) 

Defining the coefficients of the expansion in powers of g"^ as 

E{g, s) = Eo{s) + E,{s)g' + E2{s)g^ + 0{g' 

we have that 



x+(0) x-(0) 



Ms'') . 
(4.27) 

(4.28) 



Eo{s) 



dv d 


i i 


j_fj 271 dv 


y+1 ^-i 



Fo{v)-4{L-2) + 0{s 



= 41ns + 47i5-4(L-2)ln2 + o(s°) 
On the other hand the order g'^ of E{g, s) is given by 



EAs) 



dv d 
27T dv 



v+^ v-^ 



dv d 
_b 2tc dv 



+ (^-f) 



Fo(t;) + 8(L-2) + o(s°). 



(4.29) 



(4.30) 



The first line of (14.301) gives 



/"'' dv d 


i i 


F^{v) = - 


dv 


i i 


27r dv 


v-\_ 




V + - V — - 



-7T[\ns + -fE-{L- 2) \n2] [ dv ^ ^ 

j + v^ cosh^ 



nv 



^ dv 
-b^ 



2i 



2i 



L 
2 



ip [ iv ] + ip [ — \- iv 



o[s 



--7i^[ln s + - (L - 2) In 2] - 2LC(3) + o(/ 



(4.31) 
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Fo{v) + 8iL-2) 



+ 8(L - 2) + o(.°) = 6LC(3) - 14C(3) + o{s') 



On the other hand, the second hne of (14.301) gives 
dvd_ 

Adding up the two contributions we obtain 

Ei{s) = -^7r2[ln s + jE-{L-2)\n2]+ 4LC(3) - 14C(3) + o(s°) . (4.32) 

Analogously the third order in the energy, E2{s), given by 

1 1 



Eo s) 




(4.33) 



i^f (^) + 

Fo{v) - 32(L - 2) + o(s°) 



after similar and lengthy calculations, equals 

11 



Eo s) 



97r^[lns + 7s-(L-2)ln2] + 
45 



TT 



+ :-(4-L)C(3) + (62-21L)C(5) + o(s°) 



(4.34) 



Collecting all these terms we have 



2 1 1 



E{g,s) = [lns + 7£;-(i:-2)ln2] (4-^7rV + ^7rV) +/[4LC(3)-14C(3)] + 



+ g'^{4 - L)C(3) + ^7^(62 - 21L)C(5) + 0(/) + o(.°) . 
We recognize in (14.351) the expansion of the universal scaling function f{g), 

/(^7) = 4-^7rV + ^vry + 0(/), 
and remark that, consistently with fl4.2ip . relation (14.351) can be written also 

E{g, s) = ^h{g, s)f{g) + E^^'*™(^, s) + 0(/) + o{s') , 
where h{g,s) is given by (14.151) and 

E-*™(^7, s) = -2g'La3) + ^7rV^C(3) - 7rVC(3) + 9Lg\{5) + 0(/) . 



(4.35) 



(4.36) 



(4.37) 



(4.38) 
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We remark also that expansion fl4.35p agrees with result (3.16) of |16j . 

Finally, using similar techniques one computes, up to the order ^f^, the eigenvalues 
of all the conserved charges, 

1 / 1 \ ' — 1 



Qr{g,s) 



1^ 

k=l 



X+{Uk) 



1 



X (Uk) 



with r even, r > 4. Defining 

Qr{g, S) = Qr,o{s) + Qr,l{s)g^ + Qr,2(s)^7^ + 0(/) 

we obtain the following results 



Qr,o('S) 



2(-l)i-i((r-l) 



r — 1 



[(2-2^-^)L-2(l-2^-^)] + o(s 



(4.39) 



(4.40) 



Q,,i(s) = 4(-l)iC(r)[lns + 7i?-(iv-2)ln2] + 

+ L(-l)f (r + 2 - r^^)C{r + 1) + 2(-l)^(2'^+i - l)C(r + 1) + o(s°) , 

QrM = i(-l)^+^L(r + 2)(r + l)rC(r + 3)+r(-l)5+i(L-3)C(2)C(r + l) + 
o 

+ 2(-l)i-5C(r){2C(2)[lns + - (L - 2) ln2] + 7C(3) - 2LC(3)} - 
- r(r + l)(-l)5[lns + -fE - {L - 2) ln2]C(r + 2) + 

+ i(-l)3+i(r + l)(r + 2)LC(r + 3) + 

+ 2(-l)i+i(r + l)[lns + 7B- (L-2)ln2]C(r + 2) + 

+ ^(-l)^(r + 2)[-LC(r + 3) + (L - 2)(2^+3 - l)C(r + 3)] + o(.°) . 

As far as we know these expansions are new results. 

4.2 The NLIE in the s +oo hmit 

As in the one loop case, in the limit s +oo the NLIE satisfied by the counting 
function reduces to the linear equation satisfied by the density of roots, i.e. the BES 
equation. 

The counting function satisfies the NLIE 



Z{u) 



-2L arctan 2u — iL In 



- 2i 



dv 


In 1^ 




Lb 2^ 





L-2 



+ 2 arctan(M — 



h=l 




+ 2 1 —aTctan(u-v)^[Z(v)-2L(v)] 
b2n dv 



d 

dv 

z 



[Z{v) - 2L{v)] + 



2x'+{n)x-{nW) I , {^)^ 

' + i9{u,ul'} 



2x-{u)x+{u'j^'') 



(4.41) 
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where refers to the holes present inside the interval [—6,6]. 

As in the one loop case, we go to the limit s +00. In this limit, at the leading 
order in s (i.e. In s) we can drop the second and the last term in the rhs of the NLIE, 
as well as also the terms containing L{v) = Imln[l + e*^(^~*°)]. We end up with the 
linear equation 



Z(u) 



-2L arctan2'U + 2 / — arctan('U — f 

_5 27r dv 

«2 



(4.42) 



dv 


In 1^ 









1 - 



2x+(u)x~' (v) 



+ i9{u, v) 



d 



L-2 



— Z(v) + 2} arctan( 
dv 



2x~ {u)x+ (v) 

It seems natural to split the solution Z{u) as 

Z{u) = Zo{u) + Z"{u 



u — u 



h > 



h=l 



(4.43) 



where Zq{u) is the solution of the one loop part of fl4.42p . already written in (13.241) and 
Z^{u) is the solution of the higher than one loop part of (14.421) . 



z^'M 



/•^« dv 


In 1^ 









1 - 



2x+(u)x~{v) 



+ i6{u, v) 



+ 2 



^ 2^ 



2x (u)x'^(v) 

arctan(M — f ) — z In 



dv 



(4.44) 



2x+ {u)x~ (v) 



1 - 



+ 6{u,v) 



2x (u)x'^{v) , 



d 

dv 



Z^'i 



Let us now define the densities 
Id 



Zo{u) 



-po[u) 



2ns du ' s' 
In terms of them, (14.441) reads as follows 

1 



d 



2'iTsg'^ du 



Z"{u) 



(4.45) 







2ti(t^{u) 



d 



dv 



11 + [u- v) 



2i 



bo 



cT"iv) + ^ I dv4- 
du 



bo 



.d 

I— In 

du 

1 - 



2x'^{u)x {v) 



1 - 



(4.46) 



2x-{u)x+{v) 



In 



2x+ {u)x~ (v) 
2x-(u)x+{v) 



+ i6[U, v) 



This equation coincides with the BES equation. In particular, if we drop the dressing 
factor, it reduces to the ES equation - see (65) of [17]. 



4.3 An alternative derivation of the anomalous dimension 

In the large spin limit, for the leading O(lns) contribution, f{g), an elegant proportion- 
ality holds between the all-loops energy (anomalous dimension) and the Fourier zero 
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mode of the higher-than-one-loop density of roots [21]. Here, we are going to show that 
the connection extends up to the 0(s°) order. 

Regardless the normahsation, we may define the higher than one loop and one loop 
densities as a^{u) = -^F^{u), <Jo{u) = -^Fq^u), respectively, and write down their 
specific linear integral equations in the Fourier space. For we may move from fl4.1ip 
and use the fundamental Fourier transforms 



due 



-iku 



V2 



2tt h 
[±k)—e^-^Jr{V2gk) 



(4.47) 



to obtain 



sinh 



M 

2 

'OO 



+ 



sinh y J-oo \ti\ 

OO OO _ IM 

+ 2^^c,,,+i+2.(^7)(-l)"+"^(jr-i(v^(7A;)J,+2.(V2(7/^) - 

r=2 u=Q 

- Jr-i{V2gh)Jr+2u{V2gk)y^ [a"{h) + ao{h) + 27r(L - 2)] + o(s°) . 
Performing now the limit ^ 0^, we gain 



lim a^ik) = lim a"(k) = 



^MV2gh)e-^[a''ih)+aoih)+2n{L-2)]+o{s') 

(4.48) 

On the other hand, if we re-write the energy (14.121) up to the order s° by means of the 
Fourier transforms (I4.47p . we also obtain 

E{g, s) = --^ 1^ ^j^^^gh)e-'^[a^{h) + a,{h) + 27r(L - 2)] + o(s°) . (4.49) 



Upon comparing (14.481) and (14. 49 p . we size the desired relation 

a^(0) = Ve(^,s) + o(.°). 



(4.50) 



Now, in order to check this result we can repeat the perturbative expansion (I4.22p 
and extract from it the quantities Eq{s) and £"1(5). In fact, (I4.24p yields, at order (7^, 
the Fourier transform 

^i^^) = J-^[vrLfc2 + 47r|fc|(lns + 7i=;-(i^-2)ln2)] , (4.51) 

2sinh(Y) 

which just verifies 

af (0) = 7r£o(s) . (4.52) 
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Similarly, fl4.25p implies 




+ 



7T\k\^Eo{s) 

4 



+ 



6 



+ 



vr' 



k\Eo{s) 
6 



+ 



+ 27r|A;|[7C(3)-2LC(3)]}. 



(4.53) 



Finally, for k = 



af(0) 



71' 




27r[7C(3)-2LC(3)]=7ri?i(s). 



(4.54) 
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Of course, these results, (14.521 H75^ . agree with the general one (14.501) . 

5 Summary 

In this paper we have developed and applied the technique of the NLIE on intervals 
sketched in [13]. This new formalism allows to treat magnon scattering matrices with 
general dependence on the rapidities and states with roots on intervals of the real line (or 
even of complex lines). Therefore, it seems perhaps more indicated than the historical 
method presented in [Hj, if we want to study Bethe equations appearing in the context 



We have given an explicit application of the NLIE on interval for the Bethe Ansatz 
type equations describing the s/(2) sector of A/" = 4 SYM. In particular, we have written 
the exact equations which allow us to define the relevant functions (i.e. the forcing term 
F and the kernel function G) entering the NLIE. Then, we passed on to studying the 
limit of large number of Bethe roots (or spin s). In this limit and at the leading order 
Ins, the NLIE as well as the equation for the derivative of the forcing term naturally 
becomes the BES equation. If we take into account also the sub-leading correction 
0(s°), the forcing term satisfies a modification of the BES equation, with a different 
inhomogeneous term. Therefore, this equation is suitable also for a (non-perturbative 
and) strong coupling study, even if possible results in this direction should be corrected 
by eliminating the wrapping effects. Interestingly, we noticed that in the formalism 
of the NLIE on intervals the non-linear terms are negligible as s -|-oo. Therefore, 
in order to determine the eigenvalues of the conserved charges up to the order s^, it 
is sufficient to consider in their general expressions only linear terms, involving the 
forcing term. This conceptually and practically enormously simplifies their calculations 
and could suggest applications of our formalism to cases in which wrapping effects are 
under control or absent. 
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of AT = 4 SYM. 
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A Numerical evaluation of the non-linear integrals 

In this appendix, we will analyse the large spin s behaviour of the non-linear integral 
function in (12.51) and of the non-linear integral in f l2.17p . Both involve the logarithmic 
function (of Z{u)), L{u), as defined in (12.101) . In principle, we could pursue an analytic 
saddle point evaluation as in the sine-Gordon case [8], since anew, because of its forcing 
term, the counting function Z scales like the size In s, thus implying , roughly speaking, 
a correction of order ~ e~^°* = we will see that this easy conclusion should be not 
far from reality. Nevertheless, we find here more instructive for us to gain some flavour 
by performing a still easy numerical exercise. In fact, for simplicity's sake we restrict 
ourselves to the one loop case and, in the aforementioned spirit, we approximate Zo[u) 
in Lq{u) by the forcing term Fq{u) as given by ( ]3.52p . In the numerical implementation, 
we must learn to write the integrals for finite e: thanks to the Cauchy theorem, we may 
move the integration contour in the complex plane as far as we do not meet a singularity 



with (pQiu^v) defined in (13.81) and &o('S) determined by the condition (13.71) . which gives 
6o(s) ~ 0.472 s. In fact, if we should assume exactly &o(s) = s/2, we would observe an 
external hole jumping into the interval of integration at a certain value of s. Of course, 
this would give a deceiving discontinuous dependence of the integrals (lA.ip on s. In 
other words, the value s/2 is not a good one as for the separator [H] between the roots 
and the external hole. Once clarified this crucial point, we may evaluate the integrals 
with Mathematica. For instance, we report in Fig. 1 the behaviour of flA.ll) with s on 
the X-axis, for m = 1 and e = 0.1: the leading ~ 1/s decrease can be easily spotted. 
Eventually, we want to numerically evaluate the non-linear integral in the one-loop 
energy, namely the last hne of (I2.17p . once we speciahse the observable 0{v) = q2{v). 
We still need to take into account the "lateral" contributions due to the integration 
over {—ie,ie) and then, upon using the leading order solution Go{v,w) of Section 3.3, 




(A.l) 



(f)o{u,bo{s) 
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Figure 1: Comparison between numerical evaluations (dots) of 6Zq{u = 1; s) and their 
best fit as ai/s + 02/5^ + a^/s^ (line), with ai ^ 0.001, 02 ~ —0.105, 03 ~ —8.137. 



we obtain 



-Im / — q2{v) I dw— [Go{v,w - ie) 

J-bo{s) J-bo(s) dv 

6{v-w + i e)] In [l + e'^^^""-''^] + 
dy { d 



r-bo(s) 



+ 



(A.2) 



TT 



[ dv 



[Gq{v, -bo{s) - iy)- 
d 



S(v + bo(s) + iy)] ^ In h + e^^o(-^oW-*?')l 
dy 

^ [Go{v, bo{s) - zy) - 6{v - bo{s) + zy)] + e^^"(^"(«)-^) 

dv dy 



1 



Actually, since the counting function is odd, we may substitute the f-derivative of 
Go{v, w) directly the antisymmetric combination (13.371) . The dependence of SEq{s) on s 
is plotted in Fig. 2, where we still have e = 0.1, and also in this case the behaviour seems 
to be suitably fitted by a polynomial of 1/s, yielding again a 1/s leading contribution. 
Unfortunately, this one loop set-up cannot exhibit the appearance of the logarithms 
In s, in agreement with the large s expansions of [HI [201 [23] • 
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